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Abstract 

Work by Mandal and Jain [S. S. Mandal and J. K. Jain, Solid State Commun. 118, 503 (2001)] 
suggests that interaction induced mixing with the second composite fermion Landau level can 
lead to renormalization of the electron correlation function exponent in the fractional quantum 
Hall effect. In the work reported here a similar mixing with the second electronic Landau level 
is studied in the v = 1 integer case. The ground state is calculated by use of the Hartree-Fock 
algorithm, and the electron density and electron correlation function on the edge are calculated. It 
is shown that the interaction gives rise to oscillations in the density profile. In particular, a short 
range interaction gives a profile qualitatively similar to the results reported by Mandal and Jain. 
On the other hand, no renormalization of the correlation function exponent is found. 
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Introduction 

In the past few years there has been renewed interest in the edge physics of the QH 
system initiated by the discrepancy between the observed tunneling resistance [2, 3] and a 
simple application of Wen's chiral Luttinger liquid description of the edge [4, 5]. According 
to this description the Luttinger parameters will have universal values determined by the 
topological properties of the bulk. For a Luttinger liquid the electron correlation function, 
| (ip ji (x)ip(x')) |, exhibits a power law behaviour, \x — x'\~^, in the large \x — x'\ limit, where 
7 is determined by the Luttinger parameters. In the simplest case, v = l/(2m + 1), the 
correlation function exponent will have a universal value determined by the filling fraction; 
in the v — 1/3 case it can be shown that 7 = 3. The observed tunneling resistance does 
indeed confirm a Luttinger liquid description of the edge in the sense that the asymptotic 
behaviour of the edge correlation function exhibits a power law, but the exponent deviates 
from the universal value predicted by the general theory. 

Numerical works by Mandal and Jain (MJ) [1, 6] and Tsiper and Goldman [7, 8] have 
attributed the discrepancy to the effect of the electron-electron interaction. Their results 
suggest that the the exponent is not universal, but is renormalized by the interaction. How- 
ever, a different view is taken by Wan, Evers and Rezayi [9], who suggest that the discrepancy 
between experiment and theory can be attributed to the details of the edge confinement. In 
their study, they also found a change of the correlation function due to the interaction, but 
only in the case of a so-called 'hard edge', which was used by MJ and Goldman and Tsiper. 
When a 'soft edge' was introduced, or a neutralizing background potential was introduced 
in the 'hard edge' case, the universal exponent was retained. They suggested that the expo- 
nent is indeed universal, and that the edge confinement is responsible for the experimental 
results. 

Mandal and Jain [1] studied a system in the circular gauge and introduced an interaction 
between composite fermions (CF) that induced mixing with the second CF Landau level. 
In the CF description of the fractional QH effect, the ground state of the system is made 
up of CFs fully occupying one or more CF Landau levels; i.e. the fractional QH effect is an 
integer QH effect of CFs. The electronic ground state is then found by multiplying the CF 
ground state with a Jastrow factor and then projecting onto the lowest Landau level (LLL). 
MJ introduced a Coulomb interaction and found the ground state by exact diagonalization 
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in a Hilbert space spanned by a Laughlin state and states with a single CF excited to the 
corresponding angular momentum state in the second CF Landau level. This corresponds to 
a 'hard edge' since scattering to higher angular momentum states is not allowed. They found 
that the interaction gave rise to oscillations in the density profile and to renormalization of 
the asymptotic correlation function exponent. A Yukawa-type interaction also renormalized 
the exponent, but with a value different from the Coulomb case. They suggested that the 
exponent is not universal and its value is determined by the details of the interaction. 

We see no obvious reason why a similar renormalization effect should not be present in 
the v — 1 case if the electron-electron interaction causes mixing with the second electronic 
Landau level. To see whether such an effect is present, we study a model with linear 
geometry. Each single particle orbital in a completely filled LLL is allowed to mix with the 
same momentum orbital in the second Landau level (2LL). Scattering to higher momentum 
states is not allowed, which corresponds to a 'hard edge' condition, as in ref. [1]. To 
determine the the ground state wavefunction of the interacting system, we use the Hartree- 
Fock (HF) algorithm, rather than exact diagonalization. MJ's results suggests that the 
interaction form is not important for a renormalization to take place. For that reason 
a Gaussian interaction is chosen here, a choice which gives an easy control of both the 
interaction length and interaction strength. 

An advantage of working in the v — 1 regime is that a calculation involving a large 
number of particles can easily be done. This, combined with the linear gauge, means that 
large system lengths can be achieved since the system length is proportional to the number 
of particles and kp l . The calculation reported here involves 1001 particles and the Fermi 
momentum is kp = lOi^ 1 . This gives rise to a system length of about 300 magnetic lengths. 
The choice of linear geometry gives rise to two edges. Effects of the edge-edge interaction 
in the v — 1 regime will not be discussed here since it has been studied previously [10]. 
The cited work only included states in the LLL and scattering to higher momentum states 
were allowed, which makes the work complementary to the study reported here, where 
only mixing with the same momentum state in the 2LL is allowed. To mimic the positive 
background potential found in real systems a harmonic potential is introduced, and it is 
used to compensate for the expansion of the system due to the repulsive particle interaction. 

The electron density and the electron correlation function are calculated in the ground 
state. We find that that the interaction gives rise to oscillations in the density profile. For a 
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short range interaction the density, close to the edge, qualitatively resembles the edge profile 
reported in [1] and [8]. For the correlation function no renormalization of the asymptotic 
exponent is found. 

Model 

As already discussed we will work in the linear gauge where the vector potential is given 
by A x = —yB, A y = 0. The single particle Hamiltonian can then be written as 

where it is assumed that eB > and the magnetic length, is — \/h/eB, and the cyclotron 
frequency, oj c = eB/m, have been introduced. The Hamiltonian, (1), is explicitly transla- 
tionary invariant in the x-direction and energy eigenstates that are also eigenstates of p x can 
therefore be defined. With the momentum quantized as as p x = Hk, the Hamiltonian is re- 
duced to a one-dimensional harmonic oscillator with potential minimum at the /c-dependent 
position y k = —l 2 B k. The energy eigenstates can be written as 

*Pkn(x, y) = Ne lkx i) n (y - y k ), 

where Af is an appropriate normalization factor and ip n is the n'th energy eigenstate of a 
harmonic oscillator. Note that the system is degenerate with respect to k. 
Let us now introduce a harmonic background potential in the ^-direction: 

h = l-ip-eAy + l^y 

1 2 1 2 ( W c \ 1 iJ 2 .. 2 

= T^Py + ^1 V + —V. + 2 
2m y 2 y muj z c J 2mu) c 

In the last line an effective cyclotron frequency, ZJ C = eB/m = ^Juo 2 + uo 2 , with an effective 
magnetic field B = B 2 + muj 2 /e 2 , has been introduced. We see that the first two terms in 
(2) can be identified as the Hamiltonian of an electron in the effective magnetic field B. Since 
the x-momentum is still a constant of motion it can be quantized, and the energy eigenstates 
will have the same form as before, except for a modified magnetic length £b = y fr/eB and 
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with y k = — ^-£ B k. The last term in (2) lifts the degeneracy in k and the energy eigenvalues 
are given by 

e kn = Ahw c [n + - + A 



2 J V^c/ 2m ' 

where a dimensionless variable has been introduced 



a = (e B /e B ) 2 = zj c /u c = yji + (co/u c y. 

If we assume periodic boundary conditions of period L in the x-direction, the energy 
eigenstates corresponding to the LLL and the 2LL can be written as 




xexp - h/ + 4V 7^2 ), (3) 



where n = 0,1. The wavefunction (3) is centered at y — —£ 2 B k/A 2 , which means that we 
have a correspondence between the ^/-position of the particle and the particle momentum 

hk. 

Let us introduce the electron interaction with a Gaussian form, V(r) = V e~ a2r2 , which 
allows for easy control of both the interaction strength, Vo, and the interaction length, 1/a. 
We will assume the ground state to have the form; 

kp 

|GS>= n (^(^40 + ^(^)4)10), (4) 

k=— kp 

where c kn creates a particle of momentum hk in the n'th Landau level and Co(k) and C±(k) 
are mixing coefficients. Normalisation of the ground state requires that |Co(/c)| 2 + |Ci(/c)| 2 = 
1. This form easily conserves total linear momentum because each single particle orbital 
only includes states with the same /c-number. In the weak interaction limit the mixing with 
the 2LL will be small and we can assume that C (k) « 1 and %^ C 1. In an expansion 
to first order in ^4|j we see that the ground state is a superposition of a filled LLL ground 
state and states with one particle promoted to the 2LL. The restriction that all momenta 
has to be less than kp, corresponds to a 'hard edge'. 

The single particle orbitals will be of the form <^fc(r) = J2n=o l Cn(^)V ; fcn(r). From symme- 
try about the x-axis and the correspondence between k and y we would expect the orbitals 
to behave as \<fk( x , y)\ — \V-k(x, —y)\. We see from (3) that ipko(x,y) is symmetric and 



ipki(x,y) is antisymmetric (except for the phase factor e tkx ) when (k,y) — > (—k,—y). We 
would therefore expect Co(k) to be symmetric and Ci(k) to be anti-symmetric when k — > —A;. 
In terms of the ground state (4) the total energy of the system is given by 

E GS = (GS|tf|GS> 

kp j fc F 

= 5^ (Vil'iki) + « {(<Pi<Pj\V\<Pi<Pj) - (<Pi<Pj\V\(pj(pi)} , (5) 



2 

i=-k F i,j=-k F 



where 



((Pi\h\(pj) = J drip*(r)hipj(r) 



and 

{^ 3 \V\ Vm ) = J dTdr'tf(r)<p*( r >)V( r -r')Mr)<Pi(r') 

Minimizing the energy (5) leads to the Hartree-Fock (HF) equation which determines the 
single particle orbitals 

k F ~\ 

h+ £ (J, - K ^ \ ^( r ) = HQMr), (6) 

q=-kp ) 

where h is the single particle Hamiltonian and X(k) is an eigenvalue. The operator I k is 
defined as 

h = J dv' V l{v')V{v-v') Vk {v'l (7) 
and the exchange operator K k by 

K k 4>(r) = J dr' V * k (r')V(r - rV(r> fc (r). (8) 

By introducing <^fc(r) = J2 n =o i C n {k)ipkn(?) into the the HF equation and multiplying with 
/ dripl n (r) and using the fact that the V'fcn's are orthonormal, the HF equation takes the 
form 

F^k) F n (k)) ycm) \ci(k) ) 

where 

kp 

(k) 

q=—kp ni,n2=0,l 

x {{kn; qni\V\km; qn2) — {kn; qni\V\qri2, km)} (10) 
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It is easily shown that F nm (k)* = F mn (k), which has been used in the derivation of (9). The 
matrix element in (10) is defined as 

(hrii, k 2 n 2 \V\k > l n' l ; k' 2 n' 2 ) 

= J ^ ini (r)4 n2 (r')y(r-r> M< (r)^ n ,(r'). (11) 

In the limit L^> £b the x-integrals can be performed. By introducing the gaussian interac- 
tion the matrix element takes the form 

(kirii, k 2 n 2 \y\k' x ri! x \ k' 2 n 2 ) 

= ^(2A)(-^^^)/ 2 5 fcl+fc2 , M+ ^e-^-^) 2 /^ 



x / dydy'e-^y-y'^ 2 

x(y' + £ B k 2 /A 2 ) n2 e^' +eBk ^ A ^ A / 2 
x (y + £ B k[/A 2 ) n '^ +iBk '^ 2A / 2 
x (y' + £ B k' 2 /A 2 )^e-^ +eBk '^ A2 ^ 2 , 

where the integration variables y and y' are dimensionless. 

The problem of finding C(k) = (Co(k), Ci(k)) has now been reduced to a 2 x 2 eigenvalue 
problem by equation (9). The problem has to be solved iteratively: One starts by calculating 
the matrix ¥(k) from the non-interacting ground state, C(k) = (1,0). Eq. (9) is solved for 
C(k) and C(k) corresponding to the lowest eigenvalue, X(k), and highest eigenvalue, \(k), 
respectively. The new ground state is then constructed from the set of all C(k) and the 
procedure is repeated until convergence of C(k) is achieved. For the final set of all C(k) to 
constitute a ground state of the system, the corresponding energy Eqs has to be lower than 
the energy E quq2 of an excited state constructed by replacing one state C(q 2 ) with C(qi). 
The requirement is therefore that T qim = Eqs — E qim > for all q\ and q 2 , where T qim 
can be written as 

2 

r 91l?2 = A( gi )-A(g 2 )- Y, KM C MCn 3 (q2)C ni (q2) 

K}=i 

x {(gini; q 2 n 2 \V\qin^ q 2 n A ) - (qiUi, q 2 n 2 \V\q 2 n 4 ] g 2 ^4)} • 

Due to the translational invariance in the x-direction the electron correlation func- 
tion will only depend on the relative distance in the x-direction, C(x;y,rj) = 



(ip^(x : y + 77/2)^(0, y — 77/2)). In terms of the ground state coefficients C n (k) and by using 
the expression for the field operator, ip(x, y) = J2 n J2 k ipkn(x, y)ckn, the corrrelation function 
is given by 

kp 

C(x; y, V )= £ { C *o( k Mo& V + v/2) + C{(k)%l)* kl (x , y + V /2)} 

k=—kp 

x {C (k)4> ko (0, y - 77/2) + d(A;)^i(0, y - ry/2)} . (12) 
The electron density is found from (12) by setting x = rj = 0; 

p(y) = C(x = 0;y, V = 0). (13) 

Results 

For the numerics we have chosen k F = lOf^ 1 and iV = 1001 particles, which gives a 
system length of L = tt(N — l)kp X « 300£_b. From the correspondence between the position 
in y-space and fc-space we see that the width of a free system is about W = 2k pi's — 20£ B . 
To easily relate the interaction parameters to the system size, we define the interaction 
length as a — 1/a, and it is given in units of magnetic lengths. For the rest of this paper 
all length scales are given in units of the magnetic length and all wave vectors in units of 
inverse magnetic length. The interaction strength, Vo, is given in units of Jtw c , while the 
harmonic oscillator frequency is given in units of u> c . 

The ground states in terms of C (k) and Ci(k), the electron density (13) and the electron 
correlation function at the the edge, C(x;y = £ 2 B kF,i] = 0), have been calculated for three 
different cases with interaction lengths given by a = 1£b, o = 5£ B and a = 10£ B . 

To make it easier to compare results with a non-interacting system, it is preferable for 
the interacting system to have a system width of the same size as the non-interacting case. 
At the same time the requirement that T qitq2 > 0, for all q\ and 52, has to be fulfilled to 
ensure that the ground state corresponds to a true minimum. To meet these demands both 
the frequency of the harmonic background potential, u, and the interaction strength, Vo, 
have been tuned in such a way that both demands are satisfied. 

To find the number of HF- iterations required, C n {k) has been plotted, as a function of 
k, for different numbers of iterations and for a limited number of particles. The lowest 
number of iterations that does not show a change in C n (k) from one iteration to the next, 
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FIG. 1: The figure shows the density profiles for a system entirely in the LLL with uj = and 
the interacting cases with different interaction lengths. The density is measured in units of the 
non- interacting bulk density, l/2ir£ B . 

is chosen. This number of iterations have then been used for the full HF-calculation with 
1001 particles. The number of iterations for the a = 5£b and a = 10£b case is 5, and for 
the a = 1£b case 9 iterations have been used. 

Figure 1 shows plots of the electron density for a non-interacting system and the three 
interacting cases with different interaction lengths. We see that the interaction gives rise to 
oscillations that reach across the whole system, where the oscillation length increases with 
increasing interaction length. Close to the edge the a = 1£b case qualitatively resembles 
the edge profile reported in [8] and [1] for the u — 1/3 case. Also the oscillation length, 
roughly 3£b, is of the same magnitude as in the cited references. When the the interaction 
is increased beyond a = 10£b, no qualitative differences can be observed in the density 
profile. A possible explanation is that the interaction will effectively behave like a constant 
background potential, when the interaction length is comparable to the system width. 
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FIG. 2: Plots of the mixing coefficients Co and C\ as a function of k, in units of £g , for the cases 
shown in Figure 1. The number of HF-iterations is 5 for the a = 5£b and a = 101 b cases and 9 for 
the a = 11b case. As pointed out earlier we see that Co(fc) is symmetric about k = and C\{k) is 
antisymmetric. 

Figure 2 shows plots of the mixing coefficients Co and C\ as a function of k for the three 
cases. For the a = 1£b case we see that the oscillations found in the density profile are also 
present in the mixing between the Landau levels. It is also evident that the relative mixing 
between the Landau levels is of the same magnitude, even if the variation in Vq between the 
systems is quite large. 

The similarity of our a = 1£b density profile and the results reported in [8] and [1] suggests 
that we might see a similar renormalization of the exponent of the correlation function at 
the edge. Fig. 3 shows plots of the correlation functions for the different cases, and we see 
that this is not the case: The three different cases approach the same asymptotic behaviour 
as the non- interacting case. A linear fit to the log- log plots for a distance from ~ 8£b to 
~ 40£b shows agreement with the asymptotic behavior of the non- interacting system with 
7 = 1. 

Conclusions 

The study presented here is based on a similar construction as by MJ for v = 1/3, but 
for v = 1 and with a Gaussian interaction. The ground state, for three different values 
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FIG. 3: The absolute vale of the correlation function at the edge, \C(x;y = l 2 B kF,r] = 0)|, as a 
function of x for the systems in Figure 1. The correlation function is measured in units of the 
non- interacting bulk density, l/2ir£^. 

of the interaction length, was calculated using the HF-agorithm. The density profile for 
a = 1£b shows oscillation effects similar to what had been reported by MJ and by Tsiper 
and Goldman [8]. On the other hand the correlation function exponent shows no sign of a 
renormalization due to the interaction. MJ's results suggest that the interaction form is not 
important for a renormalization effect to take place. In this respect the difference in results 
between the v — 1 and v — 1/3 case is more likely to be due to fundamental differences 
between the integer and fractional systems, than to the difference in interaction type. We 
have no simple explanation for this difference. 
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